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OSCILLATORY COMBUSTION OF A LIQUID-OXYGEN JET WITH GASEOUS HYDROGEN 

by Marcus F. Heidmann 
Lewis Research Center 

SUMMARY 

Time variations in combustion parameters are measured in a zone of oscillatory 
combustion to aid in formulating a model of the oscillatory-combustion process. The 
study centers on a single 0. 028- inch-diameter liquid-oxygen jet reacting with diffusely 
injected gaseous hydrogen in a two-dimensional circular combustor. Instability is 
coupled with the first spinning transverse mode. 

A correlation of pressure and gas velocity shows a predominantly in-phase relation 
approaching that of a plane traveling wave. Gas velocities agree with calculated values 
for an isentropic wave process using experimentally determined gas properties. 

Photographs of the liquid jet correlated with phase position in the wave show the 
response of the jet to its changing environment. Mass removal from the jet is con- 
tinuous but varies with time and is in the form of unresolved small drops or vapor. It 
appears to lag behind the pressure and velocity change. Mean jet length agrees with that 
predicted for steady transverse gas velocities. 

Radiation surveys indicate energy release at and beyond the position of mass release 
from the jet. The distributed radiation pattern may be attributed to gas-phase mixing 
processes. Time variations in gas radiation are basically in phase with pressure and 
primarily correspond to an isentropic process for a gas composition in thermodynamic 
equilibrium. Cyclic reaction beyond this in-phase and concentration-dependent behavior 
is indicated by a skewing of the radiation wave shape in the direction of increasing time. 
A portion of the energy released appears to lag behind the pressure change. 

A conceptual model of instability based on the observed parameter variation is 
presented. Energy addition to the wave is dependent on the mass released from the jet, 
gas mixing, and a concentration-dependent reaction process. The effect of mass re- 
lease is amplified by a steady gas flow. 


INTRODUCTION 

An ultimate goal in combustion research is to describe analytically and to inter- 
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relate the various processes encompassing the combustion phenomena. In steady com- 
bustion, such descriptions are frequently time or space functions in which time and 
space are in themselves related. In oscillatory combustion, time and space parameters 
may vary independently, and analytical relations not only become more complex but may 
vary substantially from those for steady combustion. 

Transient and oscillatory combustion processes have been studied by various in- 
vestigators. Shock-tube experiments provide a large source of such data. Studies of 
the combustion process (refs. 1 and 2), of chemical kinetics (refs. 3 and 4), and of the 
separate process of atomization (refs. 5 to 7) are examples of shock-tube applications. 
Analytical- model studies (refs. 8 to 10), although conducted to obtain stability criteria, 
indicate the parameter variations to be expected from an interaction of processes. The 
complete combustion system has also been used to study the interaction of processes 
(refs. 11 to 13). Such an experimental approach is used in this investigation. 

The purpose of this study is to inquire into the time variations of combustion param- 
eters at various positions in a zone of oscillatory combustion, to examine these param- 
eters for deviations from steady processes, and to formulate a model for oscillatory 
combustion. The mass and energy release processes are examined in particular be- 
cause they provide the driving force for sustained oscillations. 

The investigation is made in a two-dimensional circular combustor with radial flow 
of propellants and combustion gases. An isolated 1-inch-diameter combustion zone 
centered on a 0. 028-inch-diameter oxygen jet reacting -with hydrogen is observed mac- 
roscopically for mass and energy release processes during instability. The acoustic 
environment is that of traveling transverse oscillations. A single condition of instability 
is studied. The oscillatory pressure ratio is 1. 91 at a frequency of 4350 cps for the 
first spinning transverse mode. 
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jet breakup length 

initial jet breakup length 

combustor pressure 

nitrogen orifice throat pressure 

specific gas constant 

radius 

gas temperature 
combustion gas velocity 
liquid velocity 
nitrogen gas velocity 
mean tangential velocity 
wave velocity 

combustion gas weight flow rate 
nitrogen gas weight flow rate 
ratio of specific heats 
gas density 
liquid density 


COMBUSTOR CHARACTERISTICS 

The experimental data were obtained from the two-dimensional circular combustor 
shown in figure 1. Operating characteristics of this combustor are described in refer- 
ence 14. The combustor cavity is 8 inches in diameter and 1/2 inch in length, and 
propellants and combustion gases flow radially from the circumference to an exhaust 
nozzle in the center of one cover plate. The injector consists of 48 radially alined 
oxygen orifices and indirect passages for gaseous hydrogen. Although the configuration 
differs from usual combustor designs, the orientation between wave motion and propel- 
lant flow simulates axial-flow-combustor characteristics. 

Figure 2 shows the details of the injector manifolding. Oxygen orifices are 0. 028 
inch in diameter and have a length-diameter ratio of 6. 5. Hydrogen is introduced 
through slots located on both sides of the oxygen orifice. The hydrogen path is indirect 
and spreads and diffuses upon entry into the cavity. The hydrogen slots occur only at 
the angular position of alternate oxygen orifices. This deviation from complete 
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symmetry resulted from development 
changes that were terminated when an 
unstable combustion device was obtained. 
Observations were made through the 
window ports shown in figure 1. The 


I | | I ^ I. y //////////////?< ! — region of observation is a 1-inch- 

^ diameter circle including an oxygen 
/ /^mk//////^^ I orifice and a hydrogen injection slot. 

■^Oxygen manifold Rad — I Fused quartz windows recessed 1/16 

0.488 

inch from the wall were used. 

<t 

A circular flow of combustion gases 

Figure 2. - Injector configuration. (Dimensions in inches. )_ ^ induced tangential motion of a 

steady flow of nitrogen gas. This circu- 
fnjeXn^^vx lar flow was required to establish a 

/ \ self-sustained spinning transverse mode 

/ \ of oscillation. Nitrogen gas was in- 

Exhaust/ — | jected diametrically opposite the win- 

nozzle / dow port and did not enter the field of 

\ / view because of combined radial and 

\ / circular flow, as can be seen in fig- 

X '\j/Vindow/~ > i ure 3. Normal test firings were 

l. 5 seconds in duration with nitrogen 

Figure 3. - Mean path of nitrogen flow. injected after 1 second of steady burn- 

ing. The instability was programed for about 0. 3 second. 

Combustion parameters at a single operating condition are investigated. The condi- 
tion is a mean chamber pressure of 260 pounds per square inch absolute, a metered 
oxygen- hydrogen weight flow ratio of 4. 5, a total propellant flow rate of 0. 7 pound per 
second, and a nitrogen flow rate of 0. 22 pound per second. This condition is in an in- 
stability region where pressure ratios of about 2 occurred at a frequency of about 
4400 cps. 


Figure 2. - Injector configuration. (Dimensions in inches.) 


Nitrogen^ 

injection 


Exhaust/ 
nozzle l 


Figure 3. - Mean path of nitrogen flow. 


MEASUREMENT TECHNIQUES 

Combustion parameters during instability were obtained by the use of the following 
experimental techniques. 

Chamber Pressure 


A crystal- type transducer was used to indicate chamber pressure. The transducer 
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in a water-cooled mounting is shown in 
figure 4. It was located 1/2 inch from 
the circumference - the same radial 
distance from the exhaust nozzle as the 
midpoint of the window port. The 
pressure signal was displayed on a 
cathode-ray oscilloscope and photo- 
graphed. The dynamic-pressure in- 
strumentation was supplemented by that 
used for performance evaluation. 


High-Speed Motion Pictures 


Combined silhouette and self- 
luminous photographs were taken of the 
liquid- oxygen jet with a high-speed 
motion- picture camera. Figure 1 shows the optical arrangement. Light from a concen- 
trated zirconium- arc lamp was passed through a filter so that visible wavelengths longer 
than 5800 angstroms were transmitted through the combustor. The silhouette of the 
oxygen jet formed by this red light was photographed on high-speed color film balanced 
for 3200° K light; hence, only the red emulsion of the tricolor film was used for this 
image. The characteristic radiation of hydrogen-oxygen combustion is reported in 
reference 15 to be concentrated at short wavelengths. The spectral distribution and 
film sensitivity were such that combustion radiation predominantly exposed the blue 
emulsion of the film. The green emulsion was not used. The superpositions of 
transmitted- and combustion- luminosity photographs on the same film were examined 
separately by the proper selection of filters. Oxygen-jet characteristics and the com- 
bustion zone surrounding the jet were obtained from such films. 

In the photography of oscillatory combustion the camera framing rate was adjusted 
to a value near the oscillatory frequency. This gave a stroboscopic effect in which 
successive pictures were taken at small but progressively changing phase positions in the 
pressure cycle. Exposure time for each photograph was about one-tenth of the wave 
period. 



Figure 4. - Details of high-frequency pressure transducer. (Dimensions in 
inches.) 


Streak Photographs 

Streak photographs were taken for displacement, velocity, and luminosity data. 
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The photographic arrangement was similar to that 
for motion- picture photography (fig. 1). Trans- 
mitted and emitted combustion light was recorded 
on color film with a motion-picture camera. Dis- 
placement data were recorded as viewed through 
a 1/32-inch-wide slit aperture located on the 
combustion side of the window. Slit and camera 
orientation were adjusted to indicate time varia- 
tions along either radial or tangential paths. 

Time displacement was displayed on the 
film of both the heterogeneously luminous com- 
bustion gases and the silhouetted oxygen jet. 
Streaks produced by the combustion gases were 
used for gas velocity indications in the vector 
direction of the slit orientation. Absolute evalua- 
tion of velocity, however, may not be precise 
because of the limitations in continuous tracking 
and identification of the gas zone. Streaks of the 
silhouetted oxygen jet provided additional infor- 
mation on jet behavior during instability. 

Approximate radiation intensities were also 
obtained from the streak photographs. Film transmission densities for blue light were 
obtained from recording densitometer measurements. Time variations in transmission 
densities at a specific position were obtained. Density measurements were converted 
to relative radiation intensity from a calibration of the film. The calibration curve 
obtained by time and aperture changes of exposure is shown in figure 5. 
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Figure 5. - Blue-emulsion calibration of color film. 


Time Correlation 

Transducer signals were time correlated by being simultaneously displayed on an 
oscilloscope screen. Chopper amplifier techniques for multibeam display were used. 

A high- intensity lamp with 1-microsecond light pulses at 0. 01-second intervals provided 
a signal for time correlation of oscilloscope signals with combustion photographs. The 
light was directly photographed on the combustion film records and indicated on the 
oscilloscope with a photoelectric circuit. Measurements taken at different positions in 
the combustor were time adjusted from a knowledge of the wave propagation velocity. 
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Figure 6. - Oscilloscope record of combustor pressure. 


RESULTS AND DISCUSSION 
Chamber Pressure 

An oscilloscope record of the 
chamber-pressure oscillation is shown in 
figure 6. The oscillation in pressure is 
typical of that experienced for all runs; 
however, a detailed makeup of the wave- 
form was not completely reproducible. 
Higher harmonic components were always 
evident; they were caused by either pickup 
ringing or higher modes of oscillation within the combustor. Electronic filtering was 
not used because of the indefinite phase shift caused by such techniques. 

The chopper amplifier used in the indicating circuit gave a 26-point resolution per 
cycle. Figure 7 shows the wave shape after numerical smoothing of these data points. 

A maximum pressure of 357 pounds per square inch absolute and a minimum pressure 
of 187 pounds per square inch absolute are indicated; these values give a pressure ratio 
of 1. 91 and an amplitude that is 0. 654 that of chamber pressure. The wave shape has 
the characteristics of the pronounced symmetrical spike attributed to strong spinning 
transverse waves in reference 16. 

Chamber-pressure oscillations were measured at a combustor radius of 3-^ inches, 
whereas photographic observations were made at radial positions of 3 to 4 inches. The 
oscillatory pressure amplitude is theoretically related to radius by a first-order Bessel 
function. A variation in amplitude of only 7 percent is predicted for the window port 
area. This agrees with the experimental measurements of reference 17. 

The measured oscillatory frequency in chamber pressure is 4350 cps. An expres- 
sion for the frequency of the first spinning transverse mode (ref. 18) is 

f = 0. 293 | 

A theoretical speed of sound of 4980 feet per second has been calculated for equilibrium 
combustion products for a mixture ratio of 4. 5 with nitrogen dilution. Since the 
wave velocity is added to any bulk motion of the gases, a correction for circular 



Figure 7. - Reduced combustor pressure 1/2 inch from ( 
cumference. Frequency, 4350 cps. 
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(a) Tangentially alined slit. 
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(b) Radially alined slit. 

Figure 8. - Typical streak photographs of oscillatory combustion. 


velocity produced by tangentially injected nitrogen must be made. Experimental meas- 
urements to be presented in the section Gas Velocity indicated this velocity be be 235 
feet per second at the center of the window port. Corrections for this velocity give a 
predicted theoretical frequency of 4577 cps or an experimentally indicated speed of 
sound of 4720 feet per second. The indicated speed of sound is lower than theoretical. 
Incomplete combustion near the circumference, heat-transfer losses, and the effect of 
strong oscillations (ref. 16) could account for the difference. 


Gas Velocity 

Gas velocity determinations were made at the midpoint of the window 1/2 inch from 
the injector. Typical streak photographs for both tangentially and radially alined slits 
are shown in figure 8, and the indicated velocities are presented in figure 9. Both the 
tangential and radial velocity indications are shown correlated with time and pressure. 
The general characteristic shown is that of a spinning transverse mode. The tangential 
component is in phase with pressure, and the radial component is 90° out of phase. 

The accuracy of the velocity determinations from the streak photographs is variable. 
With the original films the accuracy for peak velocities is estimated to be within 5 per- 
cent. A greater inaccuracy occurs at lower velocities and when simultaneous tangential 
and radial displacements obscure the record. 
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Positive direction of motion 




The tangential velocity is basically 
symmetrical about a mean value of 235 
feet per second with oscillations of about 
±750 feet per second. This mean value or 
constant tangential velocity is attributed 
to the circular flow induced by tangen- 
tially injected nitrogen. It is a specific 
property of this test and not a general 
characteristic of the spinning mode. 

The circular flow induced by tangen- 
tially injected nitrogen can be theoreti- 
cally approximated from a momentum 
balance: 


100 L 



.4 .6 

Time, msec 



u 0 (w r +w n ) =w n u n + (P n t -PJ A n g 


n 


n n 


n, t 


n c 


(i) 


Figure 9. - Indicated gas velocity 1/2 inch from injector corre- 
lated with chamber pressure. 


The nitrogen was injected from a 1000- 
pound-per-square-inch-absolute source 
with critical flow through a 1/8 -inch- 
diameter orifice. These values, together 
with a propellant flow of 0. 7 pound per 
second, a nitrogen flow of 0. 22 pound 
per second, and a combustor pressure 
of 260 pounds per square inch absolute can be used to compute a circular flow. A 
velocity of 394 feet per second is required to balance the momentum equation. The 
measured value of 235 feet per second is only 60 percent of this value. The loss of 
momentum is attributed to wall impingement, jet expansion angle, turbulence, and 
friction losses. Acoustic streaming, a circular flow accompanying a strong wave that 
rotates opposite to the wave (refs. 16 and 19), would also reduce the experimentally 
observed value. 

The radial velocity is similarly biased by a mean gas flow. In this case it is caused r 
by the flow of combustion products toward the nozzle. A value of 95 feet per second was 
determined from the mean slope of the particle streaks. The instantaneous values of 
radial velocities are questionable because of indistinct streaks. Maximum values of 
380 feet per second in the direction of the circumference and 280 feet per second in the 
direction of the nozzle were indicated. The mean radial velocity of 95 feet per second 
and the indicated speed of sound of 4720 feet per second are best satisfied by a model in 
which combustion is limited by oxygen vaporization, as in that of reference 20. They 
agree with a condition where, 1/2 inch from the circumference, more than 50 percent of 
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the oxygen remained to be burned. 

In a vector sum of the velocities shown in figure 9, the contribution of the radial 
component is relatively small. The wave has the predominant characteristics of a plane 
traveling wave at the position of measurement with velocity variations similar to the 
tangential component. This agrees with analytical treatments (refs. 16 and 18) where 
plane wave characteristics exist at the circumference and deviate substantially from 
this condition only on approaching the center of the cavity. 

The oscillation in velocity is related to the pressure oscillation from acoustic 
theory. An approximation of this relation can be obtained by assuming the wave to be 
plane and one -dimensional at a point near the circumference of the cavity. For such a 
wave, the energy of a gas particle with respect to the wave is 


(u c - W) 2 
2 


+ c p T = 


( 2 ) 


If small -perturbation theory applies, the relation becomes 

(W - u c )Au c = c p AT (3) 

For an isentropic process such that 


T = CjjP^VV 

c = — gR 
p y - 1 & 

a 2 = TgRT 

a relation between particle velocity and pressure can be formulated: 

2 AP 

(W - u ) Au = * £- (5) 

V P c 

The wave velocity W as approximated by the frequency of the oscillation is 
W = 2irrf - u„ and W - u = W. The perturbation in particle velocity is then related to 
combustor parameters by the relation 
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( 6 ) 


A AP c 

Au c = ?(2irrf - u 0 ) “P^" 

The experimentally determined values of these parameters as previously presented are 
a = 4720 feet per second, f = 4350 cps, and u^ = 235 feet per second. For an average 
y of 1. 2, the velocity-pressure relation at the radius of the window port is 

AP 

Au = 2410 - 

P c 

For a pressure amplitude AP /P of 0. 654 the velocity perturbation Au is 1570 feet 

L c c 

per second or ±785 feet per second about a mean value. The agreement between the 
experimental value of ±750 feet per second and the value computed for indicated gas 
properties is relatively good. The wave appears to have isentropic properties. 


Oxygen-Jet Characteristics 

A sequence of high-speed motion pictures of the liquid jet taken from the 1-inch- 
diameter window is shown in figure 10. The relative phase position with respect to the 
pressure oscillation and the approximate magnitude and direction of the velocity vector 
are shown below each photograph. These pictures were taken in successive cycles of 
the oscillation at a repetition rate that is slightly less than the period of one cycle. 

Figure 10 and other film sequences show cyclic displacement and distortion of the 
liquid jet. The characteristics observed are the cumulative effect of time-varying forces 
caused by the velocity and pressure changes shown in figure 10. In the vicinity of the 
orifice, the jet direction approximates that of the velocity vector. At the extremity of 
the jet, a correlation becomes difficult because jet displacement and distortion are the 
cumulative effects of forces over a time interval. This time interval computed from 
jet velocity measurements is about two cycles of oscillation. 

The resolution of these photographs is not sufficient in detail to show the mass 
removal process. Mass removal appears continuous but varies with time in magnitude 
and direction; it is in the form of unresolved vapor or small drops that quickly lose 
their identity. 

These jet characteristics are a pronounced change from an equivalent stable opera- 
ting condition (ref. 14). During stable combustion, the jet length extends beyond the 
1-inch window. The oscillations reduce the length to about 1/4 inch, as shown in fig- 
ure 10. A combustion length of 1/4 inch implies a chemical conversion rate of the order 
of 0. 1 pound per second per cubic inch or 500 Btu per second per cubic inch. These 
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Figure 10. - Velocity vector and phase position for sequence of jet photographs. 


values are comparable with those for large-scale combustors of efficient design. The 
instability of this model combustor is comparable with the instability of full-scale com- 
bustors in its energy release environment, even though other design parameters may 
differ significantly. 

Figure 11 shows streak photographs of the silhouetted oxygen jet by means of 
radial and tangential camera slit alinements. Liquid displacement along a radial path 
that includes the orifice is shown in figure 11(a). The average partial slope along a 
radial path indicates a jet velocity of 30 feet per second. No significant oscillation in 
slope is apparent. Injection velocity apparently is not affected by the varying pressure 
drop across the orifice. A change in the density of oxygen streaks is evident. Such 
changes are probably caused by lateral displacement of the jet during instability. 
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Injector face 



(a) Radial slit alinement. 



Time 

(b) Tangential slit alinement 1/4 inch from injector. 
Figure 11. - Streak photographs of liquid jet in silhouette. 


Streak photographs with a tangentially alined slit (fig. 11(b)) were taken 1/4 inch 
from the orifice. An abrupt periodic change in jet cross section is evident at this 
position. A correlation with the phase position of the peak combustor pressure 
(fig. 11(b)) shows the jet disappearance to lag behind the pressure peak by 10 milli- 
seconds or more. This time lag is significant only for the position of measurement and 
would differ for other positions along the jet length. A study of the high-speed motion 
pictures indicates that the position of cross-sectional change along the jet varies during 
a cycle of oscillation. At the jet extremity, the cross section is continually distorted, 
whereas no distortion occurs at the orifice. Jet breakup as a whole is therefore con- 
tinuous but time varying and becomes abrupt only with respect to specific positions 
near the orifice. 

If mass removal or jet breakup is the integral effect of an applied force, the mass 
removal would be expected to lag behind the oscillations in force. Such delays have 
been studied in shock tubes (refs. 5 and 6) for a step change in force. Complete breakup 
can occur in about 1 millisecond, although the process begins at much shorter times. 
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For the oscillating flow field of this study, the force reverses direction in less time 
than the period for complete breakup. Additional information on the dynamic response 
of the jet would be needed to compare shock-tube results with this study. 

Breakup of a liquid jet with a steady transverse flow of gas was studied in the in- 
vestigation of reference 21. The distance for the initiation of breakup and the dis- 
tance for complete breakup were found to be approximated by the relations 
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and 



These correlating parameters were derived for nonvarying flow conditions. Apply- 
ing them to this study gives relatively good agreement when integrated mean values of 
cross velocity and gas density are used. If D Q = 0. 028 inch, = 75 pounds per cubic 
foot, Pg = 0. 1 pound per cubic foot, u^ = 360 inches per second, and u c = 4000 inches 
per second, an initial breakup distance of 0. 075 inch and a complete breakup distance 
of 0. 22 inch will result. Jet appearance in figure 10 generally agrees with such lengths. 

Oscillations in jet breakup length are expected during instability, and the streak 
photographs (fig. 11) indicate they exist. The oscillation in jet length is not obvious 
in the jet photographs of figure 10; however, an extreme condition can be estimated. 

If the jet is assumed to be completely shattered by only the positive velocity peak, the 
jet can grow no longer than the jet travel distance in one cycle. In such an extreme 
case, jet breakup length would vary by less than 1/8 inch. 

A detailed analytical description of the jet characteristics is not possible from this 
study. Jet displacement and distortion vary in both degree and direction, and a complex 
expression in time and space would be required to describe them analytically. 


Gas Radiation 

Gas radiation measured at the center of the window port is shown correlated with 
time and chamber pressure in figure 12. Radiation and pressure are basically in phase; 
however, the radiation waveform is skewed in the direction of increasing time with 
respect to the pressure waveform. This difference implies a degree of nonlinear be- 
havior in the processes relating pressure and radiation. 
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An interpretation of such radiation 
measurements is aided by the studies 
reported in references 15 and 22. In 
these studies the radiation from hydrogen 
and oxygen reactions in rocket combus- 
tors was studied spectrally under a 
variety of reacting conditions. The 
spectral distribution was attributed to 
line and band spectra of the hydroxyl 
radial (OH), molecular oxygen (O 2 ), and 
water vapor (H^O) with a weak continuum, 
which remained proportional to OH radi- 
ation. Radiation intensity peaked in the 
region of stoichiometric mixture ratios 
and exceeded the radiation of the fuel- 
rich mixtures (oxygen-fuel ratios <4. 0) by an order of magnitude. Radiation emitted by 
these reacting gases was found to be in local thermodynamic equilibrium. No evidence 
of chemiluminesence or radiation from intermediate species was noted. The radiation 
from reacting zones was considered the same as that from reaction products. 

The blue-sensitive emulsion of the color film (ref. 23) used in this study was 
selectively exposed by the OH and O 2 radiation with an exposure primarily indicative 
of OH radiation. The most intensely exposed zones in a photograph of these hetero- 
geneous gases are assumed to be gases at or near stoichiometric proportions. A cor- 
rect exposure of such zones would give little or no exposure for end products at a mix- 
ture ratio of 4. 5. For the short path length through the combustor, 5/8 inch, radiation 
is assumed proportional to concentration. On the basis of these considerations, the 
radiation measurements are assumed to be those of OH concentration in gases near 
stoichiometric proportions and an index of the mass of propellant undergoing reaction. 

Radiation varied directly with pressure, as shown in figure 12. With radiation 
proportional to OH concentration, such a change would result from isothermal compres- 
sion alone. An isentropic temperature change, however, should occur. The tempera- 
ture change could affect the OH concentration in a unit volume of gas in thermodynamic 
equilibrium by changing equilibrium rates, as well as by changing gas density. Tem- 
perature can also affect radiation intensity according to Planck's law. When all these 
effects are included for a stoichiometric mixture of gases and the measured pressure 
oscillation, the computed curve shown in figure 13 is obtained. Relatively good agree- 
ment of the amplitudes of the measured and computed curves is shown. About 75 per- 
cent of the computed amplitude is due to the temperature effect in Planck's law and 
25 percent to the change in OH concentration in a unit volume. Such behavior would be 




Figure 12. - Radiation intensity 1/2 inch from injector at jet 
centerline correlated with combustor pressure. 
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computed radiation intensity 
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Figure 13. - Comparison of measured and computed radiation intensity. 
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Figure 14. - Difference between measured and computed radiation intensity 
compared with pressure oscillation. 
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expected from a reaction process in a zone 
where the mole fraction of reactants remains 
constant. 

The difference in the wave shape of the 
calculated and measured radiation may be 
used for additional examination of the com- 
bustion process. The differences in these 
two curves are shown in figure 14. A cyclic 
variation that lags behind the pressure oscil- 
lation by about 0. 03 to 0. 06 millisecond is 
obtained. This lag, if realistic, may reflect 
the lag observed in jet breakup or mass re- 
lease. 

Dissociation and recombination of OH may 
also affect radiation wave shape. The recom- 
bination rate of OH has been reported (ref. 3) 
as having the largest effect on the kinetics of hydrogen and oxygen reactions. The time 
constant for recombination is expected to be less than 0. 005 millisecond for the condi- 
tions of these tests. Its effect would be small for an oscillating condition with a time 
period of about 0. 2 millisecond. A small time constant may be reflected in a wave shape 
change for the nonsinusoidal pressure wave shapes encountered in these tests. Some of 
the wave distortion noted in the radiation may be caused by this effect. Radiation com- 
puted on the basis of frozen compositions or a large time constant would reduce the ex- 
pected oscillation in radiation by about 25 percent. It is interesting to note that only in 
the intermediate conditions between equilibrium and frozen composition can a phase 
shift or wave distortion occur. 

Radiation measurements taken at various positions along both tangential and radial 
paths through the midpoint of the window port are shown in figure 15. The oscillations 
are similar to those at the window midpoint. The wave shape is skewed to varying 
degrees, and the amplitude of the oscillation varies. Considerably more irregularities 
in wave shape occur away from the circumference and near the lateral boundaries of the 
window. This is typical of turbulent flames, where sporadic reaction defines the end of 
the reaction zone. 

The mean value of radiation intensity may also be used to examine the reaction 
process. The mean values along the tangential and radial paths are shown in figure 16. 

A peak value is shown near the midpoint of the window. If the level of radiation is 
assumed to be an index of reaction or energy release, it can be compared with the posi- 
tion of mass release from the jet. The jet disappears about 1/4 inch from the orifice, 
and liquid is evident about 1/8 inch on both sides of the jet-orifice centerline. Energy 



(b) Transverse. 

Figure 16. - Mean value at radiation intensity in trans- 
verse and radial direction through midpoint of window. 
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release peaks 1/2 inch from the orifice along a radial path and extends beyond the jet 
boundaries in the tangential direction. Based on this interpretation, energy release is 
distributed beyond the region of mass release from the jet and tends to confirm the 
previous deduction that less than one-half the oxygen is burned at the midpoint of the 
window port. 

The evidence of radiation beyond the region where liquid oxygen is evident suggests 
that energy is released by mixing of fuel-rich and oxygen-rich gases. This mixing 
process is not readily defined from the measurements of this study. An attempt to 
interpret the radiation patterns showed that the simultaneous effect of mass transport, 
time variations in an identifiable zone, and heterogenity could not be individually re- 
solved. It appears that such oscillatory data cannot be used directly but must be statis- 
tically averaged in some manner if a quantitative evaluation is to be made. This may be 
generalized to include most measurements made locally within the combustor because 
simultaneous time and space variations are inherent to the phenomena of instability. 


Conceptual Model 

A conceptual model of the unstable combustion process in the circular combustor 
can be formulated from the measurements and observations of this study. In this model, 
oscillatory conditions are assumed to be superimposed on a steady level of combustion 
in which liquid-oxygen jets are burning in an environment of gaseous hydrogen. 

The history of liquid oxygen injected into an established state of instability is 
visualized as follows. A constant weight flow rate of oxygen enters an acoustic field 
associated with a spinning transverse mode of resonance. Gas velocity oscillations are 
in phase with pressure, as described by isentropic wave equations for mean values of 
gas properties. The entering oxygen jet is bent by the oscillating aerodynamic forces 
imposed on it. The rate and degree of bending depends on the applied force and the 
inertia of the liquid, which causes jet displacement to lag the oscillation in force. 

The applied force on the jet also causes a distortion in jet cross-sectional area. 

This distortion grows with time and depends on the magnitude rather than the direction 
of the applied force. A limiting condition of distortion is obtained after some time in- 
terval, which is described as jet breakup. Jet length at breakup corresponds to the dis- 
tance traveled by an element of liquid oxygen along the jet axis in this time interval. In 
the process of breakup, mass is continually released or removed from the jet in the 
form of vapor and small drops. The amount of mass released varies with time because 
of the time-varying cross-sectional distortion accumulated within the jet. This is a 
linear process when jet length variations are small; however, a phase displacement 
with the velocity oscillation may be expected. 
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(a) No steady circular flow. 



Time 


(b) Steady circular flow. 

Figure 17. * Conceptual effect of 
steady circular flow on velocity- 
sensitive process. 


Vapor and small drops released from the jet 
undergo nearly instantaneous chemical conversion con- 
sistent with the local mixture ratio. The instantaneous 
mass flow rate of hydrogen-rich gases, which generate 
the force for the mass release rate of oxygen, deter- 
mines this local mixture ratio. The ratio is approxi- 
mately constant (ref. 7), although oxygen rich. 

Energy release, therefore, is only a fraction of that 
potentially available from the mass released from the 
oxygen jet. 

Energy release in this case is controlled by a 
velocity-sensitive process that is relatively insensi- 
tive to the directions of the velocity vector. Such a 
process reaches a peak output twice during each pres- 
sure cycle (see fig. 17(a)). In-phase energy addition 
to the wave tends to be cancelled by out -of -phase 
energy addition. A steady circular gas flow was 
established in this study by tangential injection of 
nitrogen gas. Such a steady flow modifies conditions 
(see fig. 17(b)) by increasing the in-phase component 
at the expense of the out-of-phase component. The 
circular flow thereby amplifies the in-phase energy 
addition to the wave. 


Only a portion of the energy released has been accounted for thus far. Energy con- 
tinues to be released after some displacement from the position of mass release in the 
mixing of oxygen-rich and hydrogen-rich gases. The mixing is a turbulent diffusion 
process with combustion, flow, and instability establishing a mean intensity and scale 
of turbulence. It is a steady process in that, locally, the mole fraction of reacting gases 
is relatively constant. Energy release in a unit volume of gas, however, varies with 
the density oscillations accompanying the wave. Energy addition, therefore, is in phase 
with the pressure oscillations. The kinetics are those of concentration- controlled 
Arrhenius-type reactions. 

The instability is sustained by the sum of the energy release controlled by mass 
release and mixing. The magnitude and phase of this energy addition to the wave 
balanced against the acoustic damping of the system establish a constant amplitude of 
pressure oscillations. 
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SUMMARY OF RESULTS 


The measurement of oscillatory combustion parameters in a 1-inch-diameter zone 
near the circumference of a two-dimensional combustor gave the following results: 

1. Pressure oscillations were those of strong spinning transverse oscillations at a 
frequency corresponding to the speed of sound for the mean gas properties. 

2. Gas velocity oscillations were predominately in phase with pressure oscillations 
and were related by isentropic flow equations. 

3. The liquid-oxygen jet was in a pseudostable condition with a continuous but time- 
varying removal of mass. Jet length agreed with calculated lengths for steady trans- 
verse flows. 

4. Gas radiation was basically in phase with pressure and primarily agreed with that 
for an isentropic process in a gas in thermodynamic equilibrium. A skewing of the 
radiation wave shape in the direction of increasing time implied a lag in a portion of the 
energy released during a cycle. 

5. Radiation surveys indicated that energy is released at and beyond the region of 
mass release from the liquid jet. The energy release pattern is attributed to gas mixing. 

6. In a conceptual model, energy addition to the wave is dependent on the mass 
released from the jet and gas mixing with an amplified contribution of mass release 
caused by a steady gas flow. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 28, 1965. 
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